(0) 
= 
6 
7O 

(0) 

—_ 

Ca) 

O 
jae 

—_ 
ok 
= 

O 

2) 

=) 

(= 
= 
7O 

0) 

— 

(ok 

(0) 

O 

O 
~ 


Journal of Virology 


23 


24 


Journal of Virology 


JVI Accepted Manuscript Posted Online 18 February 2015 
J. Virol. doi:10.1128/JVI.00338-15 
Copyright © 2015, American Society for Microbiology. All Rights Reserved. 


Murine coronavirus ubiquitin-like domain is 
important for papain-like protease stability and 
viral pathogenesis 


Anna M. Mielech , Xufang Deng , Yafang Chen , Eveline 
Kindler’, Dorthea L. Wheeler , Andrew D. Mesecar , Volker 
3,4 5,6 1 
Thiel , Stanley Perlman and Susan C. Baker 


j 
Department of Microbiology & Immunology, Loyola University Chicago Stritch School of Medicine, 
Maywood, IL, USA 


*Purdue University, Department of Biological Sciences, USA 


3 
Federal Institute of Virology and Immunology, Bern and Mittelhausern, Switzerland 


4Vetsuisse Faculty, University of Bern, Bern, Switzerland 


5 
Interdisciplinary Program in Immunology, University of lowa, lowa City, IA, USA 


°Department of Microbiology, University of lowa, lowa City, IA, USA 


To whom correspondence should be addressed: sbaker1@luc.edu 
Running title: ubiquitin-like domain stabilizes papain-like protease 


Word count abstract: 250, importance 122, text 5195 


0) 
= 
6 
79) 

i) 

= 

7) 

O 
jae 

SS 

(oR 
— 

O 

72) 

=) 

= 
= 
me) 

.D 

(oe 

(0) 

O 

O 
6 


Journal of Virology 


Journal of Virology 


25 


26 


27 


28 


29 


30 


31 


32 


33 


34 


35 


36 


37 


38 


39 


40 


41 


42 


43 


44 


45 


46 


Abstract 


Ubiquitin-like domains (Ubls) are now recognized as common elements adjacent to viral 
and cellular proteases; however, their function is unclear. Structural studies of the papain-like 
protease (PLP) domains of coronaviruses (CoVs) revealed an adjacent Ubl domain in Severe 
Acute Respiratory Syndrome CoV, Middle East Respiratory Syndrome CoV, and the murine 
CoV, mouse hepatitis virus (MHV). Here we tested the effect of altering the Ubl adjacent to 
PLP2 of MHV on enzyme activity, viral replication and pathogenesis. Using deletion and 
substitution approaches, we identified sites within the Ubl domain, residues 785-787 of 
nonstructural protein 3, which negatively affect protease activity, and valine residues 785 and 
787, which negatively affect deubiquitinating activity. Using reverse genetics, we engineered 
Ubl-mutant viruses and found that AM2 (V787S) and AM3 (V785S) viruses replicate efficiently 
at 37°C, but generate smaller plaques than WT virus, and AM2 is defective for replication at 
higher temperatures. To evaluate the effect of the mutation on protease activity, we purified WT 
and Ubl-mutant PLP2 and found that the proteases exhibit similar specific activities at 25°C. 
However, the thermal stability of the Ubl-mutant PLP2 was significantly reduced at 30°C 
thereby reducing the total enzymatic activity. To determine if the destabilizing mutation affects 
viral pathogenesis, we infected C57BL/6 mice with WT or AM2 and found that the mutant virus 
is highly attenuated, yet replicates sufficiently to elicit protective immunity. These studies 
revealed that modulating the Ubl domain adjacent to the PLP reduces protease stability and viral 


pathogenesis, revealing a novel approach to coronavirus attenuation. 
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Importance 

Introducing mutations into a protein or virus can have either direct or indirect effects on 
function. We asked if changes in the Ubl domain, a conserved domain adjacent to the 
coronavirus papain-like protease, altered the viral protease activity or affected viral replication or 
pathogenesis. Our studies using purified wild-type and Ubl-mutant proteases revealed that 
mutations in the viral Ubl domain destabilize and inactivate the adjacent viral protease. 
Furthermore, we show that a CoV encoding the mutant Ub! domain is unable to replicate at high 
temperature or cause lethal disease in mice. Our results identify the coronavirus Ubl domain as a 
novel modulator of viral protease stability and reveal manipulating the Ubl domain as a new 


approach for attenuating coronavirus replication and pathogenesis. 


Introduction 

Coronaviruses are emerging human pathogens. Severe Acute Respiratory Syndrome 
Coronavirus (SARS-CoV) caused the epidemic of 2002-2003, with ~10% case-fatality ratio (1). 
Middle East Respiratory Syndrome Coronavirus (MERS-CoV) is a pathogenic virus that was 
first identified in humans in 2012 (2). As of February 3, 2015 there have been 965 confirmed 
cases and 357 deaths (3). For SARS-CoV, the virus emerged from a reservoir in bats, replicated 
in an intermediate host (civet cats), and spread to humans. The epidemic strain of SARS-CoV 
evolved for efficient human-to-human spread (4-6). Public health measures of isolation of 
infected individuals led to the cessation of the epidemic in humans; however, SARS-like viruses 
remain in bat reservoirs (7-9). For MERS-CoV, dromedary camels are now suspected as the 
likely zoonotic source for transmission to humans since MERS-CoV sequences with 99% 


nucleotide identity to human MERS-CoV isolates have been detected in respiratory samples 
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from camels (10). Although there are reports of human-to-human transmission of MERS-CoV 
(11, 12), current strains seem to cause mostly lower respiratory tract disease and are not as highly 
transmissible as SARS-CoV (13). Other human coronaviruses (HCoV-229E, HCoV-OC43, 
HCoV-NL63, HCoV-HKU1) are endemic in the human population and are the causative agents 
of upper and lower respiratory tract disease and croup (14-17). To date, there are no FDA 
approved antiviral drugs or vaccines to fight human coronavirus-induced disease. In addition, 
the potential exists for coronaviruses to emerge into the human population from endemic 
reservoirs in bats or other animals (18). Identifying viral components critical for efficient 
replication and manifestation of disease will facilitate antiviral drug and vaccine development. 

Coronaviruses are RNA viruses that encode a replicase polyprotein at the 5’-end of the 
positive-strand genome. Upon virus entry, the genomic RNA is translated to produce replicase 
polyproteins (ppla and pplab), which are processed by virally-encoded proteases. Depending on 
the virus species, the ppla encodes one or two papain-like proteases (PLPs) and one 3C-like 
proteinase (3CLpro or Mpro). These proteases process the replicase polyproteins into non- 
structural proteins (nsps) that assemble with cellular membranes and facilitate virus replication 
(reviewed in (19)). Mouse hepatitis virus (MHV), the murine coronavirus used in this report, is 
commonly used as a model system to study the replication and pathogenesis of coronaviruses 
(20). 

The MHV replicase product nsp3 consists of multiple domains including two PLP 
domains (PLP1 and PLP2), two predicted ubiquitin-like domains (Ubl), an acidic region (Ac), 
ADP-ribose-1”-phosphatase (ADRP), nucleic acid-binding domain (NAB), coronavirus group 2 
marker domain (G2M), transmembrane segment (TM) and coronavirus highly conserved domain 


(Y) (reviewed in (21), shown in Fig. 1A). The UbI-1 domain has been shown to interact with 
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nucleocapsid (N) protein, which is important for virus replication (21, 22). Previous studies 
revealed that MHV PLP1 activity is required for processing the polyprotein at the nsp1/nsp2 and 
nsp2/nsp3 sites (23). Further, the catalytic activity of PLP1 is required for efficient virus 
replication (24). PLP2 was shown to recognize and process a LXGG motif and cleave the 
replicase polyprotein at the nsp3/nsp4 junction (25). The LXGG recognition site is similar to the 
RLRGG recognition site of cellular deubiquitinating enzymes. Lindner et al., 2005 was the first 
to predict that CoV PLPs could be multifunctional enzymes with both protease and 
deubiquitinating (DUB) activity (26). Indeed, further studies revealed that SARS-CoV PLpro, 
MERS-CoV PLpro, HCoV-NL63 PLP2, PEDV PLP2 and MHV PLP2 are multifunctional 
enzymes exhibiting protease, DUB and deISGylating (ability to deconjugate interferon 
stimulated gene 15 protein, ISG15 from substrates (27—34)) activities. CoV PLP activity is 
required for processing the replicase polyprotein and predicted to modulate the innate immune 
response by deubiquitination of signaling molecules activated by pattern recognition receptors 
such as RIG-I and MDAS (29, 31, 35, 36). Structural studies of the papain-like protease domain 
of SARS-CoV (37, 38) and more recently of MERS-CoV (39) and MHV (Chen et al., in 
preparation), show that these enzymes belong to the ubiquitin-specific proteins (USP) family of 
deubiquitinating enzymes and reveal the presence of a ubiquitin-like domain (Ubl) located 
upstream of the protease domain (38). However, the role of the Ub! domain in modulating the 
enzyme activity of MHV PLP2 is unknown. 

Here, we investigated the role of the Ubl domain (designated Ubl-2) adjacent to PLP2 for 
its function in MHV replication and pathogenesis. We generated proteases and viruses 
containing mutations within the Ubl domain and found that these mutations decreased PLP2 


activity and stability. Further, we found that these mutations resulted in decreased virus 
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replication and a marked attenuation of virulence. Immunization with the Ubl mutant virus 
protected mice against challenge with wild-type virus. Overall, our data demonstrate for the first 
time that manipulation of the Ubl domain adjacent to a viral protease can ameliorate viral 


pathogenicity in vivo. 


Materials and methods 


Cells. HEK293T cells were cultured in Dulbecco's Modified Eagle Medium (DMEM) with 10% 
fetal bovine serum (FBS) and 2% L-glutamine. DBT cells were cultured in Minimal Essential 
Medium with 5% FBS, 2% L-glutamine, and 10% Tryptose Phosphate Broth (TPB). BHK-R 
cells were kindly provided by Mark Denison (Vanderbilt University Medical Center) and 
maintained in DMEM media supplemented with 10% FBS, 2% L-glutamine, and 0.8mg/ml 


G418. 


Plasmids and mutagenesis. PLP2 (amino acids 1525-1911) sequence in frame with a V5 
epitope tag was codon-optimized and synthesized by Gene Script (Piscataway, NJ) (sequence 
available upon request). Codon-optimized synthetic PLP2 sequence was cloned into pCAGGS- 
MCS vector. For mutagenesis, an overlapping PCR strategy was used with primers described in 
Table 1. The introduced mutations were verified by sequencing. The nsp2/3-GFP expression 
plasmid was kindly provided Ralph Baric (University of North Carolina). The Flag-Ub plasmid 


was kindly provided by Adriano Marchese (Loyola University Chicago). 


Protease and Deubiquitinase (DUB) Activity Assays. To determine catalytic activity of the 
PLP2 constructs, 70% confluent HEK293T cells in 12-well CellBIND plates (Corning) were 


transfected using TransIT-LT1 Reagent (Mirus) according to the manufacturer’s protocol. For 
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the protease assay, the cells were transfected with 25ng nsp2/3-GFP plasmid and 300ng 
pCAGGS-PLP2-V5 expression plasmids (wild-type and various PLP2 mutants). To assess DUB 
activity the cells were transfected with 600ng Flag-Ub plasmid, and pCAGGS-PLP2-V5 
expression plasmids (wild-type and various PLP2 mutants). At 24 hours post-transfection cells 
were lysed with 300uL of lysis buffer A containing 4% SDS and 3% dithiothreitol (DTT). 
Proteins were separated by SDS-PAGE and transferred to PVDF membrane in transfer buffer 
(0.025M Tris, 0.192M glycine, 20% methanol) for 1 h at 55V at 4°C. Following this, the 
membrane was blocked using 5% dried skim milk in TBST buffer (0.9% NaCl, 10mM Tris-HCl, 
pH7.5, 0.1% Tween 20) overnight at 4°C. The membrane was incubated with polyclonal rabbit 
anti-GFP antibody (Life Technologies) at a dilution of 1:2000 for the protease assay, or mouse 
anti-flag (Sigma) at the dilution of 1:2000 for the DUB assay. The membrane was washed 3 
times for 15 minutes in TBST buffer followed by incubation with secondary donkey-anti-rabbit- 
HRP antibody at a dilution of 1:2000 (Amersham) for protease assay, or goat anti-mouse-HRP 
antibody at a dilution of 1:5000 (Amersham). Then the membrane was washed 3 times for 15 
minutes in TBST buffer. Detection was performed using Western Lighting Chemiluminescence 
Reagent Plus (PerkinElmer) and visualized using a FluoroChemE Imager (Protein Simple). To 
verify expression of the PLP2 constructs, membranes were probed with mouse anti-V5 (Life 
Technologies) antibody at the dilution 1:5000. Mouse anti-calnexin (Cell Signal) antibody at a 


dilution 1:2000 was used as a loading standard. 


Biosensor Live Cell Assay. To determine protease activity of the Ubl mutants, the previously 
described protocol was used (40). Briefly, HEK293T were transfected with 37.5ng pGlo- 


RLKGG construct and 50ng of PLP2 expression plasmids. At 18 hours post-transfection, 
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GloSensor (Promega) reagent diluted 1:50 in DMEM (10% FBS) was added. The luminescence 


was measured using a luminometer (Veritas) every hour for 5 hours. 


Generating Ubl-2 mutant viruses. To introduce mutations into MHV A59, we used a 
previously described method (40). Briefly, plasmid encoding the MHV B subclone was 
mutagenized using primers described in table 1. Plasmids encoding the complete virus genome 
were digested with restriction enzymes, gel purified, and ligated using T4 ligase at 16°C 
overnight. The ligation reaction was isopropanol precipitated and in vitro RNA transcription was 
performed using a mMESSAGE mMACHINE Kit (Ambion) according to the following 
protocol: 40.5°C for 25 min, 37.5°C for 50 min, 40.5°C for 25 min. RNA was electroporated into 
BHK-R cells, and the electroporated cells were seeded onto DBT cells. The supernatant was 
harvested 36 hours post-electroporation and plaque assay was performed as described previously 
(24). RNA was extracted from infected DBT cells at 8 hours post-infection using RNeasy Mini 
(Qiagen) and cDNA was generated using RT? First Strand Kit (Qiagen). PCR was performed 
using replicase primers (Table 1) and purified PCR product was sequenced (amino acids 747- 


848). AM2 and AM3 were plaque purified and AM2 was subjected to deep sequencing. 


Temperature Shift Experiment. DBT cells in 6-well plates were infected with 0.1 MOI of 
wild-type MHV or AM2 at 37°C. At 2, 4, or 6 hours post infection, cells were moved to 39.5°C 
until 14 hours post-infection when supernatant was harvested and virus titer was determined by 


plaque assay at 37°C as described previously (25). 


Expression and purification of wild-type and V787S mutant proteins. The PLP2 sequence 
was cloned from pCAGGS-MCS-PLP2 (described above) into LIC vector pEV-L8, which is a 


modified pET-30 plasmid. The expression of wild-type PLP2 and V787S mutant was performed 
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using Escherichia coli strain BL21 (DE3). Cultures were grown in LB medium supplemented 
with kanamycin (50 ug/ml) at 37°C until the optical density at 600nm (A600) reached 0.6. PLP2 
expression was then induced with 0.1mM isopropyl-B-D-thiogalactopyranoside (IPTG) at 25°C 
for 6h. Wild-type PLP2 was induced at 25°C for 6 h, while the V787S mutant was induced at 
18°C overnight. Cells were harvested by centrifugation at 4,690 = g for 20 min at 4°C and stored 
at -80°C until use. Cell pellets from 1 L culture were then resuspended in 40mL buffer A (25mM 
Tris pH 7.0, 500mM NaCl, 20mM imidazole, 5nM BME) supplemented with dissolved flakes of 
lysozyme and DNase, lysed through sonication and centrifuged at 28,880 x g for 30 min (4°C). 
The supernatant was filtered through a 0.45um membrane (Millipore) and loaded onto a 5ml Ni 
HiTrap HP column (GE healthcare) pre-equilibrated with buffer A. Then the column was washed 
with buffer A supplemented with 5% buffer B (25mM Tris pH 7.0, 500mM NaCl, 500mM 
imidazole, 5nM BME) until the UV was back to the baseline. The protein was eluted through a 
gradient of 5%-100% buffer B in 30 column volumes (CV). Fractions were collected and then 
pooled after enzymatic activity and purity assessment. The n-terminal (His)s-tag was then 
removed by TEV protease (His-tagged) cleavage by incubating PLP2 and TEV-protease together 
overnight at 4°C while dialyzing into buffer C (25mM Tris pH 7.0, 100mM NaCl, 10mM BME). 
Then, free His-tag and TEV protease and uncleaved PLP2 were then separated from cleaved 
PLP2 by running the sample over a Ni’ —charged HiTrap column. The flowthrough was 
collected and then concentrated using Millipore Micron concentrators to a volume of less than 2 
mL. The concentrated sample was loaded onto a Superdex-75 Hiload 26/60 column (GE 
Healthcare) pre-equilibrated with buffer D (SOmM HEPES pH 7.0, 100mM NaCl, 10mM DTT), 


and eluted at a flow rate of 2 ml/min. Fractions containing active enzyme at high purity, as 
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judged by SDS-PAGE, were pooled, concentrated and flash-frozen in 2% glycerol using liquid 


nitrogen for storage at -80°C. 


Temperature inactivation of wild-type and V787S purified protein. Wild-type PLP2 and the 
V787S mutant protein were incubated at 25°C for different time periods (0-50 min). At each time 
point, the specific activity of both enzymes was measured at 25°C using a Synergy Multi-Mode 
Microplate Reader (BioTek) with 50 uM RLRGG-AMC as the substrate and 3 uM of each 
enzyme. The assay buffer used contained 50mM HEPES, pH 7.0, 0.1 mg/ml bovine serum 
albumin (BAS) and 2mM DTT. The experiments were performed in triplicate in a final volume 
of 100 pl using the 96-well Corning Costar black microplates. Similar experiments were carried 
out when the enzymes were incubated at 30°C. To analyze the kinetic data, the ratio of the 
reaction rate at time=t to the rate at time=0 was plotted on a logarithmic scale against incubation 
time. Kinetic data of the V787S PLP2 incubated at 30°C were fitted to a first-order exponential 
decay model (Rate,/Ratep = e“'), from which the inactivation rate constant Kingc: and half-life tj/2 
were determined. The wild-type data were fit to a line since no significant temperature 


inactivation was observed. 


Thermal melting temperature (T,,) analysis using circular dichroism (CD). Thermal melting 
analyses of the wild-type PLP2 and V787S mutant was carried out with a Chirascan circular 
dichrosim (CD) spectrometer (Applied Photophysics) equipped with a temperature control 
system (Quantum Northwest Inc.) by monitoring the CD signal at 220nm while increasing the 
temperature at a step interval of 0.4°C and at a rate 0.5°C/min. Two ml of Protein samples at 
1uM in buffer with 0.1M potassium phosphate (pH 7.5) was contained in a 10 mm quartz cell 


(Starna Cells) with magnetic stirring. Thermal scans were performed in three independent 
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experiments for both wild-type and V787S mutant MHV PLP2. The melting temperatures (Tm) 


were calculated as the first derivative peak using the program SigmaPlot. 


Generation of bone-marrow derived macrophages and virus infection. Femur and tibia from 
8-week old C57BL/6 female mice were prepared and bone-marrow was isolated, and cultured in 
IMDM with 10% FCS, 1% pen-strep, 0.1% beta-mercaptoethanol and 5% MCSF obtained from 
L929-supernatants. Additional medium was added at day 3 post-isolation and cells were seeded 
at day 6 at a density of 2 * 10° cells/well. The following day, cells were infected with wild-type 
icMHV A59 and AM2 at an MOI of 1. Medium was changed after two hours and supernatant 
was harvested at various time points post-infection for titration. Supernatant was titrated on 


murine fibroblast L929 cells (in MEM containing 10% FCS and 1% pen-strep). 


Mouse studies. C57BL/6 mice were purchased from the Jackson Laboratory. The mice were 
maintained, and experiments performed at Loyola University Chicago in accordance with all 
federal and university guidelines. 4-week old C57BL/6 male mice were anesthetized with 
ketamine/xylazine prior to intracranial injection with 600 PFU of wild-type icMHV A59 or AM2 
mutant. Weight loss and survival were monitored on a daily basis. Mice were sacrificed when 
they lost 25% of their initial body weight. In other experiments, mice infected with AM2 at 4 
weeks of age, and age-matched controls were challenged with 6,000 PFU of wild-type icMHV- 
AS59 at 9 weeks post primary infection. Survival and weight loss were monitored daily. Mice 
were sacrificed when they lost 25% of the initial body weight. For histology livers were fixed in 
4% formalin and paraffin embedded. Sections were stained with hematoxylin and eosin, and 


analyzed by microscopy. 
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qRT-PCR analysis. Brains and livers of 4-week old male mice infected with wild-type icMHV 
or AM2 were harvested and homogenized at day 3 and 5 post infection and RNA was extracted 
using RNeasy Mini (Qiagen). Reverse transcription was performed using 2 pg of total RNA and 
the RT? First Strand Kit (Qiagen) according to manufacturer’s protocol. 1 pl of CDNA was used 
to set up qRT-PCR reaction according to the manufacturer’s protocol using Single Primer Assay 
for IL-6 (SABiosciences), or replicase primers (table 1). Cr values were normalized to the 


housekeeping gene (Actin). 


Results 
Identifying the ubiquitin-like domain 2 (Ubl-2) as a modulator of MHV papain-like 
protease activity 

To identify domains and residues within nsp3 that modulate protease activity, we 
generated a synthetic, codon-optimized parent construct encoding the protease (PLP2, residues 
699-1087) in-frame with a V5 epitope tag (Fig. 1A) and introduced a series of mutations into the 
predicted ubiquitin-like domain 2 (Fig. 1B). Protease activity was evaluated using an established 
trans-cleavage assay (Fig. 1C). Our goal was to identify a domain that modulated, but did not 
inactivate protease activity. As a negative control, we generated a PLP2 catalytic mutant 
(C890A) that was inactive in the trans-cleavage assay. We generated a quadruple mutant (aa 
785-788 VDVL/SSSS) and specific single mutants: V785S, D786K, and V787S. We found that 
substitution of specific residues within the ubiquitin-like domain 2 adjacent to PLP2 consistently 
affected protease activity and reduced the abundance of cleavage product detected in the trans- 
cleavage assay when compared to the activity of PLP2-WT (Fig. 1C). To determine the degree to 


which mutants’ protease activity declined we used the recently described GloSensor assay (41). 
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We found that, consistent with trans-cleavage assay results, all Ubl-2 mutants had decreased 
protease activity. However, VDVL/SSSS, V785S, and V787S showed the highest defect in 
protease activity, while D786K showed intermediate phenotype (Fig. 1D). In these experiments 
similar levels of PLP2 wild-type and mutants were expressed in each sample. These results 
implicated the Ubl-2 domain as a modulator of viral protease activity in MHV. Further, we 
determined the deubiquitinating (DUB) activity of the Ubl-2 mutants. We found that quadruple 
mutant and the single mutants V785S, and V787S had reduced DUB activity compared to wild- 
type PLP2. In contrast, D786K mutant had similar DUB activity to wild-type PLP2 (Fig. 1E). 
Therefore, we focused on the Ubl-2 V785S and V787S mutants since these Ubl mutants 
exhibited the most significant defects in the cell-based assays. 
UbI-2 mutant virus MHV-AM72 is temperature sensitive 

To determine if substitutions in the Ubl-2 domain affected virus replication, we used 
reverse genetics (42) to introduce specific mutations into the viral genome and assessed viral 
replication. Since the phenotype of the mutations in the PLP2 assay (Fig. 1) can be different 
from the phenotype of the mutant in the context of replicating virus we generated viruses 
encoding a quadruple mutation VDVL/SSSS, designated MHV-AM1, a single mutation V787S, 
designated MHV-AM2, and another single mutation V785S, designated MHV-AM3 (Fig. 2A). 
The presence of the mutations was confirmed by sequencing within the PLP2 domain. MHV- 
AMI was highly impaired, generated pinpoint size plaques, and grew to low titer (1x10° 
PFU/ml). The low titer and pinpoint plaques generated by AM1 are likely due to the quadruple 
mutation in PLP2; however, we cannot rule out the possibility that additional mutations outside 


of PLP2 contribute to this phenotype. 
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In contrast, AM2 and AM3 viruses with single amino acid substitutions within the Ubl-2 
domain replicated to titers similar to the parent infectious clone. We found that both AM2 and 
AMs3 viruses generated slightly smaller plaques than those generated by wild-type MHV AS59. 
The small plaque phenotype can be associated with temperature sensitivity, so we evaluated 
relative plaque size of wild-type and AM2 by performing plaque assays at 39.5°C and 40.5°C. 
We found that AM2 produced small plaques at 39.5°C (about 50% reduction compared to wild- 
type), but was unable to generate plaques at 40.5°C (Fig. 2B), consistent with a temperature 
sensitive defect. AM3 also produced slightly smaller than wild-type plaques but bigger than 
AM2 plaques at 39.5°C (about 20% reduction compared to wild-type) (data not shown). Because 
the temperature sensitive phenotype was more severe for AM2 and because V787 is conserved in 
other betacoronaviruses such as SARS-CoV and MERS-CoV (Fig. 1B), we focused our studies 
on AM2. We verified the engineered sequence of AM2 by deep sequencing the viral genomic 
RNA. The 2-nucleotide substitution (gtc to ccc) of V787 was confirmed, and no additional 
changes were detected (GenBank # in progress). 

To determine the replication kinetics of AM2, we infected DBT cells with isogenic wild- 
type MHV-A59 or AM2 virus at MOI of 0.1. The cells were incubated at 37°C, supernatant was 
harvested and virus titer was determined by plaque assay. We found that at 37°C, AM2 
replicated with similar kinetics to wild-type virus (Fig. 3A). In addition, no replication defect 
was observed in bone marrow derived macrophages, important primary target cells during MHV 
infection (Fig. 3B). However, if infected DBT cells were shifted to a higher temperature (39.5°C) 
there was a progressive reduction in viral titer with increased time at the higher temperature (Fig. 
3C). When the cells were incubated at 37°C for 6 hours and then shifted to the higher 


temperature, the difference between the titers of WT and AM2 was about half a log. In contrast, 
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when the infected cells were incubated at 37°C for two hours and then shifted to 39.5°C, the titer 
of AM2 was two logs lower compared to wild-type virus. Taken together, these data show that 
AM2 encoding V787S has a temperature sensitive defect in viral replication. 
Nsp3-V787S substitution impairs PLP2 thermal stability 

To investigate the temperature sensitive phenotype conferred by the V787S substitution, 
we expressed the wild-type and V787S mutant of MHV PLP2 in E. coli and then purified the 
enzymes to homogeneity using a multi-step chromatography procedure (Fig. 4A). The thermal 
stability of the wild-type and V787S purified enzymes was then determined at 25°C and at 30°C 
by measuring the specific activity of the enzymes as a function of incubation time using the 
peptide-RLRGG-AMC as a substrate. We found that when wild-type and V787S mutant PLP2s 
were incubated at 25°C, both enzymes were stable over the incubation time course and had 
similar specific activity (Fig. 4B). The activities in Fig. 4B have been normalized to 100%. In 
contrast, when wild-type and V787S mutant PLP2s were incubated at a higher temperature 
(30°C), the specific activity of V787S PLP2 decreased rapidly over time whereas the wild-type 


protein maintained full activity (Fig. 4B). These results indicate that the V787S enzyme is 


thermally unstable at 30°C with a half-life of 27.7+1.1 min. 

To determine the relative structural stability of the purified proteins over a wider 
temperature range, we performed thermal-denaturation experiments using circular dichroism 
(CD) analysis to determine the melting temperatures (Tm) of each enzyme. Consistent with 
Figure 4B, the V787S mutant exhibited diminished thermostability (Fig. 4C). The determined 
melting temperature was 6.8°C lower for the V787S mutant protein than for the wild-type PLP2, 
suggesting that the V787S mutant PLP2 is less stable and unfolds at a lower temperature than 


wild-type PLP2. Taken together, our data indicate that the V787S substitution diminished the 
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overall protein thermal stability resulting in decreased enzymatic activity at elevated 
temperatures. 
AM? exhibits reduced pathogenesis and elicits a protective response in mice 

To determine the effect of the Ubl-2 domain mutation on virus pathogenesis, we infected 
4-week old male C57B/L6 mice with 600 PFU of wild-type or AM2 virus by intracranial 
injection. We monitored weight loss over the course of infection and mice were humanely 
euthanized if they lost greater than 25% of initial body weight. We observed that while 100% of 
mice infected with wild-type virus succumbed to infection by day 7, all mice infected with AM2 
survived the infection (Fig. 5A). We confirmed that the mutation in AM2 was maintained upon 
mice infection by sequencing within the PLP2 domain of recovered AM2. 

We detected similar viral titers in the brains of wild-type or AM2-infected mice at day 3 
or day 5 post infection (Fig. 5B). In contrast we found levels of viral replicase and IL-6 were 
reduced in the livers of AM2 infected mice at day 5 post infection compared to wild-type MHV 
infected mice (Fig. 5C). There was no difference in CCL2 or IFN-a (data not shown). 
Furthermore, we analyzed livers for changes in pathology, and observed a reduction in the size 
of the lesions, hepatocyte necrosis, and inflammation in AM2 infected mice compared to wild- 
type infected (Fig. 5D). These data indicate that AM2 is attenuated compared to wild-type MHV. 

To determine if primary infection with AM2 mutant virus protected mice from challenge 
with wild-type virus, we infected AM2 immunized and age-matched naive C57BL/6 mice 
intracranially with 6000 PFU of wild-type virus at 9 weeks post primary infection. We monitored 
body weight loss over the course of infection and found that naive mice lost significantly more 
weight starting at day 2 post-infection (Fig. 5E). In contrast, AM2 immunized mice did not lose 


weight and did not exhibit other symptoms of disease upon challenge. 13-week old C57BL/6 are 
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much more resistant to infection than 4-week mice, so as expected, no wild-type or AM2 
infected mice succumbed to the infection. These results indicate that the immunization with 
AM2 generated a sufficient immune response to protect the mice from subsequent wild-type 


challenge. 


Discussion 

Modulating viral enzyme activity is a powerful approach for generating attenuated 
viruses. Coronaviruses encode a very large replicase polyprotein that contains multiple 
enzymatic activities such as proteases, ADP-ribose-1’-phosphatase (ADRP), methyl-transferase, 
exonuclease, and RNA-dependent RNA polymerase (see Figure 1). Previous studies have shown 
that inactivation of the coronavirus papain-like protease by mutagenesis of the required catalytic 
cysteine residue or inhibition with a protease inhibitor, blocks viral replication (24, 37, 43). 
These approaches demonstrate the requirement for protease activity, but represent “all or 
nothing” approaches and do not allow for the study of reduced or impaired protease activity. 
Here, we describe an approach involving mutation of the domain adjacent to the papain-like 
protease domain, termed the ubiquitin-like domain (Ubl). Structural studies revealed that this 
domain is highly conserved in coronaviruses, although the function of this domain is unknown. 
We found that mutation of the Ub! domain destabilizes the adjacent protease domain, leading to 
loss of enzymatic activity over time (half-life ~30 min at 30°C). When this mutation is 
introduced into the virus, the mutant virus replicates efficiently in cell culture at 37°C, but is 
temperature sensitive and unable to replicate at high temperature (40.5°C). Indeed, we found 


that this virus is highly attenuated in infected mice but does replicate sufficiently to induce a 
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protective immune response against wild-type virus. Overall, our studies revealed a new domain 
that can be exploited for attenuating the pathogenesis of a coronavirus. 

Ubl domains have been identified as adjacent to many viral and cellular proteases, 
particularly the ubiquitin specific proteases (USPs), a family of deubiquitinating enzymes 
(DUBs) (44). Ubl domains can be present at the C-terminal or N-terminal end of USP catalytic 
domains, or even can be inserted into the catalytic domain. In addition, USPs may contain one 
or multiple Ubl domains (45). The Ubl domain of USP14 is required for targeting this cellular 
DUB to the proteasome. At the proteasome, USP14 removes ubiquitin molecules from proteins 
that are modified with K48-linked ubiquitin, which serves as a signal for proteosomal 
degradation (46). Alternatively, the Ubl may regulate function of DUB activity by interacting 
with other cellular proteins such as GMP synthetase, which enhances the activity of USP7 (47). 

Previous studies have shown that the presence of the Ubl domain of SARS-CoV PLpro is 
necessary for efficient expression and purification of the enzyme (37). To understand the 
mechanism by which mutations in the Ub] domain of MHV PLP2 affected function, we mapped 
the mutated residues V785 and V787 into the recently solved crystal structure of MHV PLP2 
(Fig. 6) (Chen et al., in preparation). Interestingly, the overall fold of the protease and Ubl-2 
domain is conserved between SARS-CoV PLpro and MHV PLP2. Valine 787, which is the 
focus of this study, is conserved between MHV and highly pathogenic coronaviruses such as 
SARS-CoV and MERS-CoV (Fig. 1B). Valine 787 resides within a B-sheet that helps to 
maintain the structure of the Ubl. The side chains of valine 785 and valine 787 participate in the 
formation of a hydrophobic core of the Ubl, which stabilizes this domain. By mutating valine 
787 to serine, we introduced a polar side chain that is energetically unable to pack into the 


hydrophobic core because of significant desolvation penalties. Therefore, mutation of valine- 
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787 or valine 785 likely disrupts the structure of the Ubl domain thereby leading to the overall 
destabilization of PLP2 catalytic domain at elevated temperature. As a consequence, the Ubl 
mutant enzyme has decreased thermal stability that results in a loss of the specific activity of the 
enzyme over time in vitro. Furthermore, the destabilization of Ubl-2-PLP2 has an important 
impact on virus pathogenesis since we found that MHV-AM2 is attenuated in vivo. 

Interestingly, our studies revealed that the most dramatic reduction in AM2 viral titer in 
the mice was in the liver, not the brain, which was the site of inoculation. Since MHV-AS59 is 
hepatotropic, it rapidly disseminates to the liver of infected animals (19). The reduction in viral 
titer in the liver and the decreased size of lesions in the liver (Fig. 5) could be due to either 
defective viral replication or a more rapid immune response to clear the virus from the liver. 
Similar phenotype of virus attenuated in the liver but not the brain was described for MHV-A59 
ns2 mutant (48, 49). In addition, a recent study by Zhang and colleagues showed that an nsp1 
mutant virus is attenuated in the liver but replicates efficiently in the brain (50). The attenuated 
phenotype of the nsp! mutant virus is hypothesized to be due to a more robust immune response 
to the virus. 

Recent studies on MERS-CoV showed that PLpro has DUB activity and that this activity 
is important for antagonism of innate immune response in transfected cells (31, 51). In addition, 
the arterivirus equine arteritis virus (EAV) PLP DUB activity has been shown to be involved in 
inhibition of innate immune response in infected cells suggesting that PLP DUB activity might 
be important for virus pathogenesis in vivo (52). Further studies are needed to determine the role 
of DUB activity in viral replication and pathogenesis and to determine if the attenuation of AM2 


is primarily immune mediated or due to defective replication and dissemination of the virus. 
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Coronaviruses encode multiple novel enzymatic activities and previous studies support 
the concept of modifying/mutating nonstructural proteins to generate attenuated viruses and to 
identify the role of the protein in virus replication and in causing disease in mice. Mutation of the 
N-terminal replicase product nsp1 in MHV results in an attenuated infection in mice (53, 54). In 
addition, mutation of nsp! in SARS-CoV results in a decrease in antiviral signaling and virus 
titer in infected cells (55). Mutation of the catalytic residues of the ADP-ribose-1’’-phosphatase 
(ADRP) domain does not diminish virus replication in mice but reduces production of the 
cytokine IL-6, an important pro-inflammatory molecule (56-58). Reduction in inflammatory 
response upon infection with ADRP mutant may contribute to decreased liver pathology (56). 
Genetically inactivating the enzyme ExoN, which has proofreading function, results in a 
hypermutation phenotype virus that is highly attenuated, even in aged, immunocompromised 
animals (59). Mutation of the 2’-O-methyl-transferase activity in nsp16 revealed that this 
activity is required to evade viral mRNA detection by the pattern recognition receptor MDA5 
(60, 61). Further, deletion or introduction of an inactivating mutation in nonstructural protein 2 
(ns2) in MHV results in the activation of the RNase L pathway, which induces of an antiviral 
state in liver macrophages (62). All of the aforementioned studies reveal the critical roles of viral 
nonstructural proteins in virus replication and escape from immune surveillance. Our study 
reveals a role for the Ubl-2 domain in stabilization of the PLP2 catalytic domain of MHV PLP2 


and provides a new target for viral attenuation and potentially for antiviral drug development. 
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Figure legends 

Figure 1. MHV-A59 UbI-2 domain influences papain-like protease activity. A) Schematic 
diagram of MHV-A59 ORFs and the Ubl-2 and papain-like protease (PLP) domains within 
nonstructural protein 3. Expression plasmids PLP2, and the VDVL region (aa 785-788) are 
indicated. B) Alignment of the predicted ubiquitin-like domains from MHV and MERS-CoV 
with the structural information from the Ubl domain of SARS-CoV (PDB 2FE8). C) HEK293T 
cells were transfected with plasmids expressing PLP2 wild-type, PLP2 catalytic mutant CA, 
VDV UbI-2 domain mutants in the presence of plasmid expressing the nsp2/3-GFP substrate. At 
24 hours post-transfection cells were lysed and analyzed by Western blot. D) PLP2 activity in 
live-cell assay. HEK293T cells were transfected with pGlo-RLKGG and indicated PLP2 
expression plasmids. At 14 hours post infection GloSensor was added and luminesce was 
assessed hourly. E) HEK293T cells were transfected with Flag-Ub expression plasmid, and wild- 
type (PLP2-WT) or indicated Ubl-2 mutants. Cells were lysed 24 hours post-transfection and 
analyzed by Western blot. The figure shows representative data from at least two independent 


experiments. 


Figure 2. Murine coronavirus AM2 is temperature sensitive for replication. A) Summary of 
the characteristics of Ubl mutant viruses recovered by reverse genetics. B) Representative 
plaques generated by wild-type MHV (WT) or AM2 during incubation at the indicated 
temperature. DBT cells were inoculated with either WT or AM2, incubated at the indicated 
temperature and fixed and stained 48 hours post infection. Diameter of each plaque was 
measured and average diameter and standard deviation indicated for each virus; *p<0.05, 


**y<0.001. Figure shows representative data from at least two independent experiments. 
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Figure 3. Replication kinetics of AM2.. A) AM2 replication is similar to MHV-WT at 37°C. 
(A) DBT cells were infected with AM2 or WT at MOI 0.1 or (B) Bone marrow derived 
macrophages were infected with AM2 or WT at MOI 1. At indicated time points, supernatant 
was harvested and virus titer determined by plaque assay at 37°C. C) Temperature shift reduces 
yield of AM2. DBT cells were infected at the MOI 0.1 at 37°C with AM2 or WT. At the 
indicated time point post infection, cells were moved to 39.5°C until 14 hrs pi when supernatant 
was harvested and virus titer was determined by plaque assay at 37°C. The error bars represent 
standard deviation within single experiment. Figure shows representative data from at least two 


independent experiments. 


Figure 4. Ubl-2 domain is important for MHV PLP2 thermal stability. A) SDS-PAGE 
analysis of the final purified PLP2-V787S enzyme. Wild-type and PLP2-V787S mutant PLP2 
were expressed and purified from E. coli as described in Materials and Methods. B) 
Temperature-dependent inactivation of wild-type and PLP2-V787S mutant. The specific activity 
of wild-type (WT) and the V787S mutant PLP2 enzymes were measured after incubation at 25°C 
and 30°C for different time periods. The specific activities were then normalized to the activity 
at 0 min (Rate;,/Rateo: rate at time t over initial rate). Kinetic data for the V787S mutant 


incubated at 30°C were fit to a first-order exponential decay model (Rate,/Ratey = e“). The 


calculated inactivation rate constant Kina: for the V787S mutant at 30°C is 0.025 + 0.001 min! 


and the half-life (t12) is 27.7 + 1.1 min. C) CD melting curves of WT and PLP2-V787S mutant 
PLP2. The thermal stability of WT and V787S mutant PLP2 was determined by measuring the 
CD signal at 220 nm as a function of temperature at a step interval of 0.4°C and at a rate 
0.5°C/min. Three independent experiments were performed for both WT PLP2 (grey) and 


V787S mutant (black). 
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Figure 5. MHV AM72 is attenuated and generated protective immunity in mice. A) C57BL/6 
mice were infected with 600 pfu WT icMHV or AM2 mutant intracranially and monitored for 
survival over time (N=7 for each group). B) Virus titers in brains of mice infected with 600 pfu 
WT icMHV or AM2 were determined at indicated time points p.i. by plaque assay (N=4 to 5). 
Error bars represent standard deviation. C) RNA levels for viral replicase gene and cellular IL-6 
were determined using qRT-PCR in livers at day 5 post infection (N=6). RNA levels were 
normalized to actin and amount of transcript in WT infected mice was set to 1. Error bars 
represent SEM. Statistical analysis performed using Student’s ¢ test. D) Representative samples 
of hematoxylin and eosin staining of formalin-fixed liver samples at day 5 post infection from 
mice infected with WT icMHV or AM2 (N=6 for each group, magnification X10). E) C57BL/6 
mice immunized with AM2 mutant and naive age-mached controls were challenged with 6000 
pfu icMHV intracranially 9 weeks post primary infection. The mice were monitored for body 


weight loss. Error bars represent SEM.****, p<0.0001, two-way ANOVA. 


Figure 6. X-ray structures of the MHV and SARS-CoV Ubl-2 domains. A superposition of 
the Ubl-2 domains from MHV PLP2 (blue) and SARS-CoV PLpro (orange, PDB 2FE8) is 
shown as a stereoview and resulted in RMSD of 0.74 A. Side chains of MHV UbI-2 residues 


V785 and V787, and their corresponding residues in SARS-CoV are shown as sticks. 
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Figure 1. MHV-A59 UbI-2 domain influences papain-like protease activity. A) Schematic diagram of MHV-A59 
ORFs and the Ubl-2 and papain-like protease (PLP) domains within nonstructural protein 3. Expression plasmids 
PLP2, and the VDVL region (aa 785-788) are indicated. B) Alignment of the predicted ubiquitin-like domains from 
MHV and MERS-CoV with the structural information from the Ub!l domain of SARS-CoV (PDB 2FE8). C) HEK293T 
cells were transfected with plasmids expressing PLP2 wild-type, PLP2 catalytic mutant CA, VDV Ubl-2 domain 
mutants in the presence of plasmid expressing the nsp2/3-GFP substrate. At 24 hours post-transfection cells were 
lysed and analyzed by Western blot. D) PLP2 activity in live-cell assay. HEK293T cells were transfected with pGlo- 
RLKGG and indicated PLP2 expression plasmids. At 14 hours post infection GloSensor was added and luminesce 
was assessed hourly. E) HEK293T cells were transfected with Flag-Ub expression plasmid, and wild-type (PLP2- 
WT) or indicated Ubl-2 mutants. Cells were lysed 24 hours post-transfection and analyzed by Western blot. The 
figure shows representative data from at least two independent experiments. 
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Figure 2. Murine coronavirus AM2 is temperature sensitive for replication. A) Summary of 
the characteristics of Ubl mutant viruses recovered by reverse genetics. B) Representative 
plaques generated by wild-type MHV (WT) or AM2 during incubation at the indicated 
temperature. DBT cells were inoculated with either WT or AM2, incubated at the indicated 
temperature and fixed and stained 48 hours post infection. Diameter of each plaque was 
measured and average diameter and standard deviation indicated for each virus; *p<0.05, 
**y<0.001. Figure shows representative data from at least two independent experiments. 
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Figure 3. Replication kinetics of AM2. A) AM2 replication is similar 
to MHV-WT at 37°C. (A) DBT cells were infected with AM2 or WT at 
MOI 0.1 or (B) Bone marrow derived macrophages were infected 
with AM2 or WT at MOI 1. At indicated time points, supernatant was 
harvested and virus titer determined by plaque assay at 37°C. C) 
Temperature shift reduces yield of AM2. DBT cells were infected at 
the MOI 0.1 at 37°C with AM2 or WT. At the indicated time point 
post infection, cells were moved to 39.5°C until 14 hrs pi when 
supernatant was harvested and virus titer was determined by plaque 
assay at 37°C. The error bars represent standard deviation within 
single experiment. Figure shows representative data from at least 
two independent experiments. 
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Figure 4. UbI-2 domain is important for MHV PLP2 thermal stability. A) SDS-PAGE analysis of 
the final purified PLP2-V787S enzyme. Wild-type and PLP2-V787S mutant PLP2 were expressed 
and purified from E. coli as described in Materials and Methods. B) Temperature-dependent 
inactivation of wild-type and PLP2-V787S mutant. The specific activity of wild-type (WT) and the 
V787S mutant PLP2 enzymes were measured after incubation at 25°C and 30°C for different time 
periods. The specific activities were then normalized to the activity at O min (Rate/Rate,: rate at 
time t over initial rate). Kinetic data for the V787S mutant incubated at 30°C were fit to a first-order 
exponential decay model (Rate,/Rate, = e*'). The calculated inactivation rate constant kina for the 
V787S mutant at 30°C is 0.025 + 0.001 min” and the half-life (t,,) is 27.7 + 1.1 min. C) CD melting 
curves of WT and PLP2-V787S mutant PLP2. The thermal stability of WT and V787S mutant PLP2 
was determined by measuring the CD signal at 220 nm as a function of temperature at a step 
interval of 0.4°C and at a rate 0.5°C/min. Three independent experiments were performed for both 
WT PLP2 (grey) and V787S mutant (black). 
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Figure 5. MHV AWN2 is attenuated and generated protective immunity in mice. A) C57BL/6 mice 
were infected with 600 pfu WT icMHV or AM2 mutant intracranially and monitored for survival over 
time (N=7 for each group). B) Virus titers in brains of mice infected with 600 pfu WT icMHV or AM2 
were determined at indicated time points p.i. by plaque assay (N=4 to 5). Error bars represent 
standard deviation. C) RNA levels for viral replicase gene and cellular IL-6 were determined using 
qRT-PCR in livers at day 5 post infection (N=6). RNA levels were normalized to actin and amount of 
transcript in WT infected mice was set to 1. Error bars represent SEM. Statistical analysis performed 
using Student's t test. D) Representative samples of hematoxylin and eosin staining of formalin-fixed 
liver samples at day 5 post infection from mice infected with WT icMHV or AM2 (N=6 for each group, 
magnification X10). E) C57BL/6 mice immunized with AM2 mutant and naive age-mached controls 
were challenged with 6000 pfu icMHV intracranially 9 weeks post primary infection. The mice were 
monitored for body weight loss. Error bars represent SEM.****, p<0.0001, two-way ANOVA. 
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Figure 6. X-ray structures of the MHV and SARS-CoV UbI-2 domains. A superposition of the 
UbI-2 domains from MHV PLP2 (blue) and SARS-CoV PLpro (orange, PDB 2FE8) is shown as a 
stereoview and resulted in RMSD of 0.74 A. Side chains of MHV UbI-2 residues V785 and V787, 
and their corresponding residues in SARS-CoV are shown as sticks. 
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Table 1. Primers 


PLP2-CA 
SSSS 
V785S 
D786K 
V787S 
AM1 
AM2 
AM3 


Replicase 


Forward Primer 
5'G TCA AAC AAC AAC GCC TAC ATC AAC G 


5' GGC CAA TAA GAG CAG TAG TAG TTG CAC CGT CGA CGG 
5' GGC CAA TAA G AGT GAT GTG CTG TG CAC CGT CGA CG 
5' GGC CAA TAA G GTC AAA GTG CTG TG CAC CGT CGA CG 

5' GGC CAA TAA G GTC GAT TCG CTG TG CAC CGT CGA CG 


Reverse Primer 
5' CGT TGA TGT AGG CGT TGT TGT TTG AC 


5' CCG TCG ACG GTG CAA CTA CTA CTG CTC TTA TTG GCC 
5' CGT CGA CGG TGC ACA GCA CAT CAC TCT TAT TGG CC 
5' CGT CGA CGG TGC ACA GCA CTT TGA CCT TAT TGG CC 
5' CGT CGA CGG TGC ACA GCG AAT CGA CCT TAT TGG CC 


5’ CTC GCT AAT AAG AGT TCG TCC TCG TGT ACT GTT GAT GG 8’ CCA TCA ACA GTA CAC GAG GAC GAA CTC TTA TTA GCG AG 


5’ CTC GCT AAT AAG GTT GAT TCC TTG TGT ACT GTT GAT GG 5’ CCA TCA ACA GTA CAC AAG GAA TCA ACC TTA TTA GCG AG 


5’ CTC GCT AAT AAG TCT GAT GTC TTG TGT ACT GTT GAT GG 5’ CCA TCA ACA GTA CAC AAG ACA TCA GAC TTA TTA GCG AG 


5’ GAA TGC CCG GGC GGG ATT TTT GTA TCC 


5’ CCA CAA GAT CTG CCT CGG ACA AAT C 


